Finite range effects of nuclear force in intermediate energy heavy ion 

collisions 
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Within the framework of an isospin-dependent quantum molecular dynamics model, 
the zero-range 2-body part of the Skyrme interaction is replaced by a finite-range 
Gaussian 2-body interaction. From the transverse momentum analysis in the reaction 
' of system 93 Nb + 93 Nb at energy of 400 MeV/nucleon and impact parameter 6=3 fm, 

\ it is shown that the finite-range nuclear force enhances the transverse momentum of 

the reaction system and it can partly replace the momentum dependent part of the 
nucleon-nucleon interaction. 



PACS number(s): 25.70.-z, 25.75.Ld, 21.30. Fe 



I. INTRODUCTION 



The heavy-ion transport theory is one of important research frontiers in the theoretical nuclear 
physics in the last decade. Since 1980s, the beam energy has gone up to intermediate and high 
energy regions and the high temperature and high dense nuclear matter could be formed in the heavy- 
ion collisions (HIC's), through which one could investigate the nuclear matter properties at extreme 
conditions. Correspondingly, theoretical nuclear physicists have established some microscopic heavy- 
""^5 ' ion transport theory, among which the Bolztmann-Uehling-Uhlenbeck (BUU) equation M , Boltzmann 

Langevin (BL) equation 0, Quantum Molecular Dynamics (QMD) model p), Fermionic Molecular 
Dynamics (FMD) model and Antisymmetric Molecular Dynamics (AMD) model ||, have been 
used extensively and successfully. Though these models have respective advantages at different aspects, 
they have a common character, namely, the zero-range Skyrme-like effective interaction being used 
in all these models. It is well known, however, that the nuclear force is attractive at long range and 
repulsive at short range. Hence, the finite-range forces may be as more "realistic" than zero-range 
Skyrme-like forces. With the zero-range forces, some momentum-dependent terms are quadratic or 
constant leading to irrealistic behavior while those terms vanish with finite-range forces even though 
relative momenta go high ||. At usual relative momenta, below 2kp , the behavior of zero-range forces 
are still acceptable but it is always interesting to study the features of HIC's with finite-range nuclear 
forces 0. 

Based on an isospin-dependent QMD model P-[l2[ which includes the symmetry potential, isospin 
dependent nucleon-nucleon (JV-JV) collisions, Coulomb potential for protons, and isospin dependent 
Pauli blocking effects, the zero-range 2-body Skyrme interaction is replaced by a finite-range Gaussian 
interaction. As the first investigation of the finite-range effects, we study the collective flow in reaction 
93 Nb + 93 Nb at energy of 400 MeV/nucleon and impact parameter &=3 fm since the finite range of 
nuclear force represents momentum-dependence of nuclear force in the momentum space [[l3| and the 
collective flow is sensitive to the momentum dependent interaction. The calculated results indicate that 
the finite-range 2-body interaction actually enhances the transverse momentum, which implies that 
finite-range nuclear force could replace partly the momentum-dependent part of the N-N interaction 
and therefore by using the finite-range N-N interaction the momentum dependence of the nuclear 
forces is considered naturally. 
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II. MODEL AND METHOD 



In the QMD model, nucleon i is represented by a Gaussian form of wave function: 

\D-(r t) - 1 -\r-n(t)] 2 /(4L) ipi-T/H /i\ 

Performing Wigner transformation for Eq. (1), one can get the nucleon's Wigner density distribution 
in phase space: 



where the and p, represent the mean position and momentum of the ith nucleon, respectively, the 
L is the so-called Gaussian wave packet width (here L=2.0 fm 2 ). In the QMD model, then the total 
interaction potential can read 

U tot = u (2) + u (3)^ (3) 

and the two-body and three-body potentials, and can be given, respectively, by 

U(2) = E U if = E / /i(r ! P 1 ^i(>-,>-')/i(r',p',t)rfr(Jpdr' ( |p' 1 (4) 



»'(^3) ' 



and 



^ (3) = E ^ 



= E / / ^ (r ! P,^)/, (r^p^^)^(r^p'',^)^ fe (r,r',r'')drdp C ^r'dp'dr'' C ^p''. (5) 

The two-body interaction includes the local two-body Skyrme interaction V^° c , the Yukawa 
(surface) interaction VjJ uk , the symmetry energy part V^ ym , and the Coulomb interaction V^ "'. 
Correspondingly, the total two-body potential in the QMD model is as follows, 

c/( 2 ) = uffi + u Yuk + u sym + u Coul = J2 ( u l° c + u lj uk + u ij m + u Z oul )> ( 6 ) 

for the forms of U Yuk , U sym , and J/ 00 "', one j s re fereed to Refs. |J|,|l|-|l§. Refs. @]l6) give a de- 
tailed description for the Isospin-QMD (IQMD) model. In the present isospin-dependent QMD model, 
however, Pauli blocking of neutron and proton is treated respectively, namely, the Pauli blocking is 
isospin dependent. The two-body Skyrme interaction is given as following form 

Vl° c = - r,). (7) 

From Eqs. (4) and (7), one can get the two-body Skyrme potential 

"ft - E ^ oc) - E ^i^^-^ir^ (8) 

Similarly, from the three-body Skyrme interaction 

7A 3 fe ) =t 2 ^(r i -r J -)5(r i -r fe ), (9) 
one can get the three-body Skyrme potential 
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^3°) = E u Sl c) 



E 



If one defines the interaction density. 



2 33/ 2 (2^L)3 P[ 6L J ' (W) 



P\nM) = J^Y.^ Jj ^P L l (11) 



then Eqs. (8) and (10) can be approximately combined into 



C7 ioc = a( A^u )+(3 ™ )7) (12) 
Po PO 



where the po is the normal nuclear matter density, 0.16 fm 



For the momentum dependent interaction V^ ID1 , we make use of the real part of the optical 



potential parametrized in Ref. |L7|] as follows 

Vjf DI = Sln 2 [e( Pl - p. t ) 2 + l]S{ ri - r,-). (13) 

The parameters of Eqs. (12) and (13) are given in Table I, from which one can see two kinds of 
equations of state (EOS) are commonly used. One is the so-called hard EOS (H, HM) with an 
incompressibility of K=380 MeV, and the other is the soft EOS (S, SM) with an incompressibility of 
/■T=200 MeV ||. The M refers to the inclusion of the momentum dependent interaction. 
If the 2-body Skyrme interaction, Eq. (7), is replaced by a finite-range Gaussian, i.e., 

where the A represents the finite-range parameter of N-N interaction, then one can get the Gaussian 
2-body effective potential in the QMD model 

Ugr = Y U\f au) = Y h * , exp[- {Ti ~ T ^ ]. (15) 

By comparing Eq. (15) with Eq. (8), one can find that the Gaussian 2-body effective potential 

A 2 

can be reached only through replacing the L by L + in the 2-body Skyrme potential, namely, 
Eq. (8). For the zero-range Skyrme two-body interaction, Eq. (7), the short range repulsive core 
is not considered. Physically, the Eq. (14) is equivalent to simulate the short range repulsive core 
and therefore A should have the same order as radius of the repulsive core, namely, A « 0.4 ~ 0.6 
fm. The long range attractive part of the 2-body JV-JV interaction mainly comes from the Yukawa 
interaction. In order to observe the effects of A on nuclear force, we display in Fig. 1 the finite- 
range Gaussian 2-body force as a function of the distance between the two nucleons for A = 0.0, 0.3, 
0.5 and 0.7 fm. Here we set t\ — 356 MeV-fm 3 and it should be noticed that the case of A = 0.0 
corresponds exactly to the case of zero-range Skyrme two-body force since Eq. (7) is the limitation of 
Eq. (14) as A — » 0. It is indicated in Fig. 1 that the 2-body nuclear force is attractive and becomes 
weak gradually with increment of A, which imply the finite-range enhances the repulsive effect of 
nuclear force. This phenomenon is easy to understand since the finite range of nuclear force means 
momentum-dependence of nuclear force in the momentum space. 

Similarly, the zero-range 3-body Skyrme interaction, Eq. (9), can be replaced by finite-range Gaus- 
sian 3-body interaction, i.e., 

C = ^^)^[- (r '" rj)2 A + / r '- rt)2 l. («) 
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after tedious algebra, one can get the finite-range Gaussian 3-body effective potential in the QMD 
model, i.e., 

jjGau _ TT (Gau) _ ST^ ^2 

W&Hjtk i*i&hj*k 3W[2n,J(L + %)(L+ ^)]3 

x , (A 2 + 2L)[(rj - Yj) 2 + (rj - r fc ) 2 ] + 2L(rj - r fc ) 2 

Obviously, this expression is very complicated and has not the simplicity of Eq. (15). In fact, 
neglecting the A 4 and its higher order terms, Eq. (17) can be also simplified as following form, 

Gau _ T j{Gau) _ V"^ *2 



jjGau _ ST TT( Gau ) — ST 



■+ l3K ^ 33/2 \ML + 4^)] 3 

(r f - rj) 2 + (rj - r fc ) 2 + (rj - r fe ) 
6(L+f) 



xexp[- ^ .A A,/ 1- ( 18 ) 



Surprisingly, by comparing Eq. (18) with Eq. (10), one can also find that the Gaussian 3-body 
effective potential can be reached only through replacing simply the L by L + 4j- m the 3-body 
Skyrme potential, namely, Eq. (10). Generally, the zero-range 3-body part of Skyrme interaction can 
be regarded as a good approximation for the 3-body part of N-N effective interaction. Therefore, we 
only insert the finite-range 2-body part into the QMD code in the present work. 

A 2 A 2 

It should be noted that the replacing of the L by L + =|- or L + =|- is only for the two-body or three- 
body effective potential in the QMD model and this replacing dose not affect the other components of 
the QMD model. Therefore, it is not simply to modify the width of the Gaussian wave packet in the 
QMD model. The change of the two-body or three-body effective potential only comes from the using 
of the two-body or three-body finite-range Gaussian interaction (Eq. (14) or (16) ) in stead of the 
two-body or three-body zero-range Skyrme interaction (Eq. (7) or (9) ). In fact, the present method 
is very similar with that of using the well-known finite-range Gogny interaction where the finite range 
is given by the superposition of two Gaussian function with different ranges and spin-isospin mixtures. 



III. RESULTS AND DISCUSSIONS 



Within the framework of the above modified QMD model, the transverse momentum is calculated 
for the reaction of 93 Nb + 93 Nb at energy of 400 MeV/nucleon and impact parameter 6=3 fm for five 
different potential parameter sets, namely, hard EOS (H), soft EOS (S), soft EOS with momentum 
dependent interaction (SM), soft EOS with A = 0.1 (S 0.1), and soft EOS with A = 0.5 (S 0.5). Fig. 
2 gives the time evolution of transverse momentum per nucleon for the above five cases. From Fig. 2 
one can see clearly that at the initial stage of the collisions the transverse momentum is negative for all 
cases except for the case with SM, which is easy to understand since for the cases of H, S, S 0.1, and S 
0.5, at the initial stage of collisions the attractive nuclear mean field is dominant and consequently the 
transverse momentum is negative. However, for the case of SM, the initial large relative momentum 
leads to a large repulsive momentum dependent potential which balances the attractive part of nuclear 
mean field. 

With time evolution, the projectile and target begin to collide, compress, and expand, and the 
positive transverse momentum is generated. It is shown in Fig. 2 that the transverse momenta 
have saturated after 40 fm/c for all five cases. In the saturated region, one can see that the largest 
transverse momentum is observed in the case of H and the smallest in the case of S. The transverse 
momentum in the case of S 0.1 is little larger than that in the case of S. Comparing to the case of S 
0.1, the case of S 0.5 has a larger transverse momentum which is still smaller than that in the case of 
SM by an amount of 8 MeV/c. These features indicate that the potential parameter set of soft EOS 
with A = 0.5 could replace partly that of soft EOS with the momentum dependent interaction, which 
implies that the finite-range Gaussian two-body nuclear force could replace partly the momentum 
dependent part of the N-N interaction. 
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In order to compare the present calculated results with the experimental data, we display in Fig. 3 
the rapidity distribution of the transverse momentum for above five potential parameter sets. The big 
open circles included in Fig. 3 represent the experimental data which are from the Ref. |l8) . The small 
open circles linked by lines represent the calculated results. Fig. 3 (a) shows the calculated results 
for the case of H, which indicates that the calculated results are in agreement with the experimental 
data. Fig. 3 (b) corresponds to the case of SM and it is shown that the calculated results are in 
better agreement with the experimental data. In Fig. 3 (c) the calculated values linked by solid line 
and dashed line correspond, respectively, to the cases of S and S 0.1, which are not in good agreement 
with the experimental data. In the case of S 0.5, Fig. 3 (d) indicates that the calculated results 
are basically in agreement with the experimental data. These features imply that all the calculated 
results with potential parameter sets H, SM, and S 0.5 are in agreement with the experimental data 
while those with potential sets S and S 0.1 fail to reproduce the experimental data. Therefore, the 
finite-range parameter A = 0.5 fm is very reasonable to model the finite-range effect of nuclear force. 

Fig. 4 displays the time evolution of transverse momentum per nucleon for free nucleons (A=l, 
solid circles in Fig. 4) and light fragments (A=2^4, open circles in Fig. 4) for four different cases, 
namely, Figs. 4 (a), (b), (c), and (d) correspond, respectively, to the cases of H, S, SM, and S 0.5. 
In this paper, we construct clusters in terms of the so-called coalescence model Jl9| with coalescence 
parameters Rq = 3.5 fm and Po — 260 MeV/c. One should note that the mass dependence of the 
transverse momentum shown in Fig. 4 demonstrates the well known increase in magnitude for heavier 
fragments ]2p|-|22]| . This phenomenon may be because most nucleons are emitted by the hard stochastic 
collisions and hence the effect of the mean field is largely erased in the nucleon flow. This argument 
suggests that the flow of composite fragment carries more direct information of the nuclear EOS than 
the nucleon flow. From Fig. 4 one can also find that the finite-range Gaussian two-body nuclear force 
enhances the transverse momenta of fragments and could replace partly the momentum dependent 
part of the N-N interaction. 

IV. CONCLUSIONS 

A finite-range Gaussian two-body interaction is applied in stead of the zero-range two-body part of 
the Skyrme interaction in an isospin-dependent QMD model to simulate the transverse momentum 
in the reaction of system 93 Nb + 93 Nb at energy of 400 MeV/nucleon and impact parameter 6=3 fm. 
The calculated results show that the finite-range nuclear force enhances the transverse momentum of 
the reaction system and it can partly replace the momentum dependent part of the JV-JV interaction. 
The calculated results with potential parameter set of soft EOS with finite-range parameter A = 0.5 
are basically in agreement with the experimental data. Meanwhile, it is shown that the transverse 
momentum of light fragments is greater than that of free nucleon, which agrees with the results of 
many experiments. 

This work is a primary investigation of the finite-range effects of nuclear force and further consider- 
ation is worth performing. For example, the two-body forces adopt the finite-range Gogny forces and 
meanwhile the finite-range three-body force should be included. In addition, it is very interesting to 
explore the finite-range effects of other physical phenomena. 
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TABLE CAPTIONS 



Table I The parameter sets of Eqs. (2) and (3). The S and H refer to the soft and hard equations 
of state, the M refers to the inclusion of momentum dependent interaction, and the K refers to 
the incompressibility. 





K(MeV) 


a(MeV) /3(MeV) 7 (MeV) 5(MeV) 




s 


200 


-356 


303 


1.17 






SM 


200 


-390 


320 


1.14 


1.57 


500 


H 


380 


-124 


71 


2.00 






HM 


380 


-130 


59 


2.09 


1.57 


500 



FIGURE CAPTIONS 

FIG. 1 The finite-range Gaussian 2-body nuclear force as a function of the distance between the two 
nuclcons for A = 0.0, 0.3, 0.5 and 0.7 fm. 

FIG. 2 Time evolution of transverse momentum per nucleon for five different potential parameter 
sets, namely, hard EOS (H), soft EOS (S), soft EOS with momentum dependent interaction 
(SM), soft EOS with A = 0.1 (S 0.1), and soft EOS with A = 0.5 (S 0.5) for the reaction 93 Nb 
+ 93 Nb at energy of 400 MeV/nucleon and impact parameter &=3 fm. 

FIG. 3 Rapidity distribution of the transverse momentum for five different potential parameter sets, 
namely, hard EOS (H) (a), soft EOS with momentum dependent interaction (SM) (b), soft EOS 
(S) and soft EOS with A = 0.1 (S 0.1, the small open circles linked by dashed line) (c), and soft 
EOS with A = 0.5 (S 0.5) (d). The larger open circles are the experimental data and the small 
open circles are the calculated results. 

FIG. 4 Time evolution of transverse momentum per nucleon for free nucleons (A=l, solid circles) 
and light fragments (A=2~4, open circles) for four different cases, namely, hard EOS (H) (a), 
soft EOS (S) (b), soft EOS with the momentum dependent interaction (SM) (c), and soft EOS 
with A = 0.5 (S 0.5) (d). The lines included only to guide the eye. 
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